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Abstract. Weak lensing effects are known to introduce 
non- linear couplings in the CMB temperature maps. In 
inflationary scenario, the primary CMB anisotropies are 
expected to form a 2D Gaussian map, for which, the prob- 
ability distribution function of the ellipticity defined from 
the local temperature curvature matrix has a very spe- 
cific shape. I show that lenses alter significantly the shape 
of this PDF, inducing an excess of elongated structures. 
The precise functional form is computed for both the field 
points and the temperature extrema. 

These analytical results are confirmed by numerical ex- 
periments on 10x10 square degree maps. These numerical 
results allow to investigate the effects of smoothing and 
to estimate the cosmic variance. For the best resolution 
and sky coverage of the Planck mission the signal to noise 
ratio for the statistical indicators presented here is about 
3 to 6 depending on the cosmological models. A marginal 
detection should therefore be possible. 

Key words: Cosmology: Dark Matter, Large-Scale Struc- 
tures, Gravitational Lensing, Cosmic Microwave Back- 
ground 



1. Introduction 

The detection by the COBE/DMR experiment (Smoot et 
al. 1993) of the anisotropies of the CMB temperature at 
very large angular scales, above 7 deg, has open an exit- 
ing new mean of investigation for cosmology. Many ex- 
periments that are now under development will provide 
us in the near future with precious data at smaller an- 
gular scale, down to a fraction of a degree. At such an- 
gular scales the dominant mechanism that generates the 
CMB anisotropies is the primary, linear order, coupling 
of the various cosmic fluids (photons, baryons and differ- 
ent species of possible dark matter) before and during re- 
combinaison. Whether the primary anisotropies originate 
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from quantum fluctuations in an inflationary scenario, or 
from other mechanisms such as topological defects gen- 
erated in a phase transition epoch is yet unclear. This is 
undoubtlessly a major scientific goal for the coming ex- 
periments. Indeed, one of the clearest signature of the 
topological models is that they are expected to induce 
non-Gaussian primary temperature fluctuations (Pen et 
al. 1994, Turok 1996, Barnes & Turok 1996). But even in 
the case of inflationary scenario, it is possible that sec- 
ondary effects, non-linear couplings of the radiation field 
with matter, induce non-Gaussian features. This is par- 
ticularly important at very small angular scale, down to 
the arcmin scale. I am more particularly interested here 
in the static lens effects. The effects of the gravitational 
distortion on the power spectrum of the primary temper- 
ature maps have been the subject of many investigations 
in the last decade (Blanchard & Schneider 1987, Kashlin- 
sky 1988, Cole & Efstathiou 1989, Sasaki 1989, Tomita & 
Watanabe 1989, Linder 1990, Cayon, Martinez-Gonzalez 
& Sanz 1993a, b, Fukushige, Makino & Ebisuzaki 1994, 
Seljak 1996). It is now clear that the lens effect on the 
power spectrum is small, although it is probably worth 
to take it into account in a very detailed analysis of such 
measurements. In principle, for an inflationary scenario, a 
detailed analysis of the CMB power spectrum should allow 
to constrain very accurately the cosmological parameters 
(see for instance Jungman et al. 1996). But it would be 
anyway extremely interesting to be able to have a positive 
detection of the lens effect since it can potentially be used 
to constrain the amplitude of the cosmic density power 
spectrum, independently of the constraints obtained from 
the CMB power spectrum itself. It could therefore be a 
precious test for the cosmological model(s) favored by the 
shape of the CMB power spectrum. 

In a previous paper, the four point functions induced 
by the lens effects has been investigated (Bernardeau 
1997). This is the most direct quantity that can be cal- 
culated from the coupling terms introduced by the lenses. 
One thus obtains specific properties of CMB maps induced 
by the lenses. Another possible way of investigation is to 
look for cross-correlation between the CMB temperature 
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gradients and the displacement field induced by the pro- 
jected large-scale structures (Suginohara et al. 1997). The 
signal however depends on the bias properties of the galax- 
ies. It cannot be directly interpreted in terms of cosmo- 
logical parameters. 

If the four point function is indeed sensitive to the 
intrinsic depth of the lens potential wells, it is not neces- 
sarily the best indicator in terms of signal to noise ratio 
for a realistic experiment. The aim of this paper is there- 
fore to investigate new means of detecting these effects. In 
section 2 I present numerical results showing peculiar ex- 
amples of the lens effects on temperature maps. The visual 
impression is that the lens effects induce a change in the 
topological properties of the temperature maps, and not 
so much in the local temperature distribution function. 

This is the motivation for the investigation of other 
statistical indicators that can reveal the lens effects with 
a better efficiency. Particularly interesting can be the sta- 
tistical properties of the local curvature of the temperature 
maps. The matrix of the second order temperature deriva- 
tives has indeed specific statistical properties in case of a 
2D Gaussian field. These properties arc affected when the 
lens effects arc taken into account because of the induced 
mode couplings. More specifically a quantity that can be 
a good tracer of the lens effect is the probability distri- 
bution function of the local ellipticity since lenses tend 
to systematically stretch the local temperature patches. 
Compared to the 4-point function the motivation for con- 
sidering such a quantity is then double. First of all, the 
cosmic variance is expected to be smaller: It is given by 
statistical properties related to the second order deriva- 
tive, and therefore more sensitive to the small scale tem- 
perature fluctuations. Secondly the effects of lenses on the 
four-point function were found to be proportional to the 
cosine of the angle joining the observational directions (see 
Bernardeau 1997), which tends to substantially reduce the 
lens signal when it is averaged for the computation of the 
four point moment. Such a cancellation is not expected 
for the distortion effect on the curvature. The derivation 
of the lens effect for the distribution function of the local 
ellipticity is presented in Section 3. In Section 4, results 
are confronted with numerical experiments. In particular 
I estimate the cosmic variance for the ellipticity statistics. 
In the last section I discuss the dependence of the lens 
effects on the cosmological parameters. 

2. The physical mechanisms 

First of all let me recall the basic effect induced by the 
lenses. I am interested here in the static lenses caused by 
the presence along the line-of-sights of large-scale concen- 
trations of matter. These concentrations of matter are of 
course intrinsic to any theory of large-scale structure for- 
mation. The effects presented in this paper are therefore 
robust consequences of the growth of structures. The am- 
plitude of the effects, however, depends not only on the 



amplitude of the local fluctuations but also on the global 
cosmological parameters: the lens effects are all the more 
large that the CMB plane is distant from the lenses. The 
dependence of the signal on the cosmological parameters 
is discussed in the last section. The static lenses (contrary 
to the case of moving lenses, see Birkinshaw & Gull 1983) 
do not change the background temperature of the CMB, 
but only induce a displacement field so that the observed 
temperature fluctuation in a given direction 7 is the one 
coming from a slightly shifted direction, 7 + £(7), 

5 ob, (7)=(5 prim. (7 + e(7)); (1) 

where £ (7) is the displacement field induced by the lenses 
along the line of sight in the direction 7. Written in terms 
of the gravitational potential ^, the displacement reads, 

£(7) = -2/ d X — T7- r W>(7,x), (2) 

Jo ^(Xcmb) 

where \ an d 2? are respectively the radial and angular 
distances. These quantities are identical in case of a back- 
ground universe with zero curvature. 2?cmb is the angular 
distance to the last scattering surface. The lens population 
identifies with the large-scale structures of the Universe 
that are present along the line of sights. It is important 
to have in mind that the displacements induced by these 
matter concentrations are at most of the order of 1 arcmin 
for cluster cores, whereas the instrumental resolution for 
the future satellite missions 1 is at best 5 arcmin. As a re- 
sult, it is safe to neglect the influence of the critical regions 
and the lens effect is assumed here to be fully in the weak 
lensing regime. 

A consequence of this remark is that the effects in- 
vestigated in this paper can be entirely described by the 
two-point correlation properties of the lens population. 
The precise reason is that the displacement field being a 
small perturbation of the primordial field, its leading con- 
tribution can be expressed perturbatively solely with the 
two-point correlation function of either the displacement 
field or its derivatives such as the local gravitational con- 
vergence. This situation was already encountered in the 
computation of the induced temperature four-point func- 
tion (Bernardeau 1997). Therefore, the lens effects will 
depend quantitatively only on the shape and magnitude 
of the cosmic density power spectrum, P(k), defined from 
the Fourier transform of the local density field, 

/d 3 fc 
— ^ S k exp(ife • x), (3) 

with then, 

(6h6 k >) = 5n iia , c (k + k')P(k). (4) 

1 Note however that ground based interferometer instru- 
ments can potentially reach a better resolution. 
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Note that in this analysis, P(k), is the actual power spec- 
trum, that includes possible nonlinear evolution of the 
density field. Written in terms of the modes 8k the dis- 
placement field reads, 



«7) = I d X W ( X ) J ^-^ x 
i fcj_ 



(5) 



k 2 V( X ) 



6{k) exp [iD(x) fc± • 7 + ik r x] 



with, 



Ax) = 



3»p X>(xcmb-x)P(x) 
a(x) £>(xcmb) 



(6) 



The function w(x) gives the lens efficiency function for 
sources located on the last scattering surface. The second 
moment of the displacement field then reads 



rXCMB pXCMB r Ail 

(e) = J o d x w( X ) J d X ' w( X ') J 



d 3 k 

r 



kl 



P ( fc ) fci^2 ex P [mx)-V( X ') fc ± - 7 + ikr (X-X')] • (7) 

Applying the small angle approximation, the integral in 
k r introduces a Dirac function in X = x' so that 



/ d 2 fc /'(/••) 



CXCMB f J2l 



) 2 k 2 V 2 



(8) 



It is then simple to see that the displacement field can be 
expressed in terms of a 2d projected potential, p 2 d, defined 
by 



P2d(0 = I d X -^P[l/V( X )} 
with 

^ ^ [ d 2 l p 2 d(0 



(9) 



(10) 



(2tt) 2 I 2 ' 
The lens effects are thus entirely determined by p 2 d- 

2.1. Numerical Experiments 

The effect of the mechanism described in (1) can be easily 
visualized in a simple numerical experiment. A realistic 
temperature map can be generated from a given power 
spectrum, C;. A 2D displacement field can also be gen- 
erated as a realization of a Gaussian process with the 
power spectrum p 2 d obtained from a given 3D power spec- 
trum. In the following I use a standard CDM model for 
both the temperature map and the displacement field. The 
normalization of the temperature power spectrum is arbi- 
trary; the one of the displacement field will be discussed 
in terms of erg, i.e. the r.m.s. of the density fluctuations 
at 8ft -1 Mpc scale. A technical simplification is to assume 
that the displacement field is Gaussian. This is obviously 



not exact, specially at small scale, but as it has already 
been mentioned, the lens effects depend only on the power 
spectrum of the projected lens density. To mimic the ef- 
fects of the lens on the CMB maps it is thus not necessary 
to have a full 3D realization of the large-scale structures. 
Note that the non-linear effects for the evolution of the 
lens density field are implicitly taken into account with 
the use of the non-linear power spectrum P(k) (Jain et al. 
1995, Peacock & Dodds 1996, following the approach of 
Hamilton et al. 1991). The effects of the lenses are then 
computed from the equation (1): it means that the lensed 
temperature map is given by the primordial one, but on 
a now irregular grid. The lensed temperature map is ob- 
tained by a "regridding" of the temperature map, with 
a local linear interpolation of the temperature with the 
nearest neighbors. This algorithm is usually available in 
standard mathematical packages. 

The two maps of Fig 1 show the effects of lenses. The 
effect is of course magnified so that it can be easily de- 
tected by a visual inspection of a very small map. The 
value of cr 8 would be here of about 2. The size of the maps 
is 4x4 degrees, for a pixel size of 1.17'. On the lensed map, 
one can see that the temperature patches are still present 
with the same temperature contrasts. They are however 
displaced and deformed. The most striking feature is prob- 
ably that there are very sharp edges and extended area 
of constant gradients. All these features suggest that the 
local curvature is probably a good indicator of the lens 
effects. 

3. The lens effects on the local temperature cur- 
vature 

3.1. The local curvature for the primordial temperature 
map 

The temperature curvature is defined from the matrix of 
the second order derivatives of the local temperature, 



(11) 



dxidxj 



that can be written, 

d 2 S T d 2 S T 




v da:ida:2 



T + 91 92 

92 r-gi 



■ (12) 



The previous equation defines the scalar field r and the 
vector field g. Note that a rotation of the coordinate sys- 
tem with an angle 8 transforms the components of in 



K 2e ■ g, 



(13) 



where H29 is the rotation matrix of angle 29. Thus the 
vector g behaves rather like a pseudo-vector. For a Gaus- 
sian temperature field it is easy to show that r, g\ and 
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Fig. 1. An example of primordial temperature map (left) deformed by the lens effect (right). It corresponds to a CDM model. 
The figures are 4x4 square degree large. The lensed effects have been magnified to make them visible. The displacement field 
has been multiplied by about 4. 



gi are three Gaussian independent fields. At a given point 
their variances are given by, 



(r 2 ) = *h 

(al) = {al) = 1*1 

with (72 given, in the small angle approximation, by 



(14) 
(15) 



(16) 



where C; is the temperature power spectrum. This entirely 
defines the statistical properties of the quantities associ- 
ated with the primordial temperature field. 
A quantity of interest is then the ellipticity, 



where g is naturally given by, 



9 = M9i+9l- 



(17) 



(18) 



The local ellipticity expresses the fact the temperature 
fluctuations may be locally isotropic or not. Let me be a 
bit more precise. First of all it is always possible to assume 
that 52 equals zero (by a proper choice of coordinates) so 
that g = \g\\. Consider then the curvature matrix cy and 



its 2 eigenvalues. They identify with T—g and r+g. If they 
have the same sign, as it is the case for a local extremum, 
it is also the sign of r and consequently the one of the 
ellipticity. In such a case, when the ellipticity is positive, 
r — g is positive and the ellipticity is smaller than 1/2. In 
general the absolute value of extremum ellipticities is less 
than 1/2, i.e. elongated peaks correspond to ellipticities 
close to 1/2 and round structures close to whereas values 
larger than 1/2 correspond to regions of eigenvalues of 
different sign such as saddle points 2 . 

The distribution function of the local ellipticity can be 
easily computed from the known statistical properties of 
r and g. One finds that (eg. Bond & Efstathiou 1987), 



P 



prim. 



(e)de 



I de 



(l + 8e 2 ) 3 / 



2 ' 



(19) 



As there are no statistical differences at all between the 
hot and the cold regions in the primordial map, and since 
the lens effects do preserve this property, in the following 
I will consider only the cases of positive ellipticities. It 
is important to note that the distribution (19) does not 
depend on the shape of the temperature power spectrum. 
It is specific of a 2D Gaussian field. 



2 Note that although this discussion is focussed on points of 
vanishing temperature gradient, the curvature tensor and the 
local ellipticity are defined on every field point. 
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In (19) e runs from to oo (when limited to a posi- 
tive ellipticity) . For practical reason I make the change of 
variable, 

£ =l/(^ + l) (20) 

so that e runs from to 1. The distribution of e can be 
easily computed from the one in e, 

P r m (e)de= (1 _ 2 2 £ £ + d ; e2)3/2 ; 0<e<l. (21) 

This distribution peaks at e = 0.5. The numerical results 
will be discussed in terms of e because a binning of con- 
stant width is obviously more appropriate for e compared 
to e. 

The effect of lenses is, from the visual inspection of 
figure 1, to extend the regions with very asymmetric cur- 
vature, that is to favor regions of ellipticity close to 1/2. 
The aim of the coming sections is to quantify analytically 
this effect. 

3.2. The lens effects 

3.2.1. The lens effects on the curvature matrix 

We have now to consider the local curvature matrix in the 
presence of lenses. The first derivative of eq. (1) leads to 

(Sr-h = (S P T im -)A5fi+^i) (22) 

where Sfj is the Kronecker symbol. This expression makes 
naturally intervene the lens amplification matrix, Aij , 

^-"♦M-O^ 1 i-™ )• < 23 > 

where Id is the identity matrix and ip is the projected po- 
tential of the matter density fluctuations along the line-of- 
sight (the power spectrum of which obviously is P2d(l)/l 2 )- 
Note that in the weak lensing regime the matrix A is al- 
ways regular, so that the extrema are conserved by the 
lens effects (at least as long as filtering is not taken into 
account). 

The second order derivatives of eq. (1) lead to 

(S p T rim -),k^. (24) 

Two terms appear. The second term vanishes for a peak, 
and in such a case the lens effect can be written in a rather 
compact way, 

c obs. = A . C prim. . ^ ( 25 j 

where C is the curvature matrix. This expression is familiar 
in studies of weak lensing effects on background galaxies. 
It is also interesting to note that a peak cannot be trans- 
formed in a saddle point, and inversely: from relation (25) 
it appears that the determinants of C prim - and C obs - have 
the same sign, as long as the lens effect is not critical. 
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3.2.2. The local ellipticity distribution 

The computation of the statistical properties of the matrix 
element of the curvature matrix of the lensed map then 
requires the knowledge of the statistical properties of the 
lens contributions. First of all the displacement field is 
statistically independent of the temperature fluctuations 
of the CMB. Secondly the displacement is small, so that 
it is safe to do a perturbative calculation with respect to 
the displacement. It implies that only the second moment 
of the displacement, and of its derivatives, will appear in 
the final results. 

The components of the amplification matrix A can be 
written in terms of the local gravitational convergence k 
and distortion 7, 

W^"" 71 ). (26) 

V -72 1 - K - 71 J 

Similarly to the case of the curvature matrix the r.m.s. 
values of n, 71 and 72 are related to each other through, 

<7i 2 > - <7!) = \(* 2 ) - \«l (27) 

although k and ji may not be Gaussian distributed. This 
is a simple consequence of geometrical averages, that take 
advantage of the assumed statistical isotropy of the dis- 
placement field. The second moment of the local con- 
vergence, a 2 , is related to the projected power spectrum 
with 3 , 

The relation between a K and the normalization of the lo- 
cal density spectrum, a%, makes intervene the projection 
effects that depend on the cosmological parameters. That 
will be discussed in more detail in the last section. 

The second order derivatives are all independent of the 
first order derivatives. Only 4 different terms are gener- 
ated, $1,12, $1,22 and $2,22- The others are identical 
to one of those 4 because the displacement is potential. 
Then we have, 



(£2,22) 
(Ci 2 12 ) 
(?i 2 22 ) 


= 5/16 s 2 


= 5/16 s 2 


= 1/16 s 2 


= 1/16 s 2 


(£l,ll£l,22) 


= 1/16 s 2 


(£1,1262,22) 


= 1/16 s 2 



(29) 



and the other cross-correlation terms are zero. This is 
again a simple consequence of geometrical averages. The 
quantity s is also related to the projected power spectrum, 

s 2 = / -^/ 4 P2d(0- (30) 

3 This expression ignores the effects of filtering and is thus 
not realistic. I reconsider this at the end of this section. 
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The calculation of the statistical properties of the cur- 
vature elements is then now quite straightforward. One 
has to (perturbativcly) invert the relation between the 
quantities associated with the lensed map and the primor- 
dial map, and then to compute the distribution function 
r and gi from their distribution function in the unlensed 
case. These calculations are in practice a bit lengthy but 
can be performed without much difficulty with a formal 
calculator. One then has, 

P° hs - (r, 9i, 92) = P prim '(r, 91, 92) [l + Qi o* K + Q2 s 2 ] , (31) 
with 

Qi = -24 + 24 t 2 - 2 t 4 + 39 g 2 - 17t 2 g 2 - 8g 4 , (32) 
Q 2 = (-3 + r 2 + 2 ff 2 )/8. (33) 

From these expressions it is easy to compute the new dis- 
tribution function for the ellipticity, 



C(e)=PT1e) 



l + a 



2 18(-l + 20e 2 - 16 e 4 ) 
(l + 8e 2 ) 2 



(34) 



Remarkably the term in s 2 disappears. The relation (34) is 
one of the central results of this paper. In the following this 
result is compared with results of numerical experiments. 

3.3. The extrema ellipticity distribution 
3.3.1. For a Gaussian field 

The calculation of the ellipticity of the extrema has been 
examined in detail by Bond & Efstathiou (1987) in case of 
a 2D Gaussian field. The principle of the calculation relies 
on the computations of the number of peaks, N pea k(r, g) 
of a given curvature (within the range dr and d 2 g). If the 
extrema have the directions 7 P in a given sample then, 

^V P oak(T,g)drd 2 g = 

f d 2 7^J Dirac (7- 7p )Pr.(r, 9 |V(5 T -0)dTd 2 9 , (35) 
J v 

where Pr.(r, <?|V(5t = 0) is the probability distribution 
function of the local curvature under the constraint that 
the local gradient vanishes. In case of a Gaussian field this 
expression simplifies since the local gradients are statisti- 
cally independent of the curvature, i.e., 



Pr. prim '(r, 9 |V(5P rim -=0) = Pr. prim '(r, g). 



(36) 



curvature. Indeed, for peaks we can Taylor expand the 
local temperature gradient in, 



(S T )4l) = (*r)ji(7p)(7 - 7„)j 



(37) 



since, by definition, (#T),i(7p) = 0. As a result we have 
for the primordial field, 

m ( \a m J(r,g) Pr.(r,g)dT d 2 g 
A r pcak(r,g)dTd g = N tot . -j^-r, rrr— . - , 9 , (38) 



/ J(T, S )Pr.(r, g )drd2 g ' 



where iVtot. is the total number of peaks, and J(r,g) is 
the Jacobian of the transform between the positions and 
the gradients, i.e., 



J(r,g) = \Det( Cij )\ = 
As a result we have, 



\r 2 -9 2 \ 



(39) 



\r 2 -g 2 \p pr[m -(T,g) dr d 2 g, 



prim./ \ j j2„ tf l J \'itfj ^ y 5 /Ar\\ 

C k ^) drdg= /|T2 _ g2|pprim . (Tig) dr d2g , (40) 



which, after simple calculations, implies that, 

im . f )d 24V3(l-4e 2 )ed e 
P ^ lJ (l + 8e 2 ) 5 / 2 ' 



(41) 



when the integration is limited to the maxima (extrema 
with positive curvature), that is when t 2 > g 2 and r > 0. 

3.3.2. The peak ellipticity distribution with lens effects 

In this section I investigate the effects of lenses for the el- 
lipticity distribution of the extrema. The principle of the 
calculation is very similar derivation to the Gaussian field 
case. The only change is that in principle the local curva- 
ture may not be independent of the gradient. Actually for 
peaks we know that (25) is exact so that we can express 
the curvature of the lensed map in terms of the curvature 
of the primordial map for which the statistical proper- 
ties are completely defined. The effect of lenses will then 
show up in the distribution function of the local curvature 
and in the expression of the Jacobian. Taking these two 
changes into account we have, 

P&(^-») =PleZ(r,9i,92) [1 + Qzo-l] , (42) 



with 



= -28 + 24r 2 + 39.g 2 - 2t 4 - 17 t 



2 „2 



8g 4 



(43) 



Note that in this case it is necessary to make sure that 
the integral is correctly normalized. The number density 
of peaks remain unchanged with (sub-critical) lenses but 
the variances for the curvature or the gradient are pertur- 
bativcly affected. From the previous result we have, 



The trick to complete the calculation is then to express 
the Dirac delta functions in 7 — 7 in terms of the local ^P k 



Pbs. (e) 



prim. / \ 

Ppk (e) 



1 + 



2(-15 + 512e 2 - 112e 4 ) 
(l + 8e 2 ) 



.(44) 



A straightforward examination of (34) and (44) show that 
the effect of lens distortion is significantly larger for the 
peaks compared to the field points. 



Weak Lensing on CMB Maps 



7 



3.4- The effects of filtering 

So far, in these analysis, I have ignored the effects of filter- 
ing. Actually the CMB temperature maps will be observed 
with a finite resolution. The measured temperature map 
is thus given by 

~St(j) = J d 2 7 ' W{\j - 7 '|/0o) <5r(V). ( 45 ) 

where W is a window function that in practice can be 
assumed to be Gaussian. It is anyway possible to further 
filter the resulting maps. In particular it can be very inter- 
esting to filter the large angular scale modes out in order 
to reduce the cosmic variance. In all cases, the filtered 
lensed temperature map is given by 

^ bs '(7) = /dV^(l7-7'l/^)^ rim -(7' + e(7'))- (46) 

It is obviously not equivalent to the lens distortion ef- 
fects applied to the filtered primordial temperature map! 
The two operations do not commute in general. This is 
a standard problem when one is considering the effects 
of non-linear couplings. In this particular problem there is 
however a limiting case for which the problem can be com- 
pletely solved. Indeed the CMB temperature maps has an 
intrinsic small-scale cutoff due to the Silk damping effect. 
If then the filtering scale is much smaller than this intrin- 
sic scale one can assume that locally S^ 1 " 1 '^') w <5^ rim '(7) 
when f' is close enough to 7. As a result we have 

(^T bs 0,*(7)« (^ rim ')j x 

[^ + /dVW(| 7 -Y|/(9o)€(7 , ) J -,i]- (47) 

The filtering effect applies only to the displacement field. 
The same situation occurs for the 2nd order derivatives. 
In this limiting case, the effect of filtering is then simply 
to reduce the lens effect by filtering the displacement field. 
It implies in particular that the a\ factor that intervenes 
in the expressions (34) and (44) is actually given by, 

°* K = \J (0P2d(OW? o (O- (48) 

In practice however this approximation may not be very 
accurate. Even at 4.5' arcmin resolution scale (the best 
channels of Planck surveyor) about 20 to 25% of the pri- 
mordial signal is filtered out because of the finite reso- 
lution. One problem that then should be addressed with 
numerical experiment is the validity of this approximation. 

4. Numerical results 

The numerical results are based on a series of 20 realiza- 
tions of both the CMB temperature fluctuations and the 
displacement maps. Each map is initially 10x10 deg 2 large, 
with a pixel size of 1.17' (these are 512x512 grid maps). 
After rcgridding to take into account the lens effects and 



filtering, the actual useful size is 9x9 square degrees. In 
Figs. 2-4 I give the averaged ellipticity distributions ob- 
tained from these 20 realizations, and the error-bars cor- 
respond to the estimated cosmic variance of one map. 



4-1. The ellipticity distribution 



2.0 




0.0 0.2 0.4 0.6 0.8 1.0 

eps. 



4 '- 




-0.2 r 



-0.4 L_L n , , , , 3 

0.00 0.10 0.20 0.30 0.40 0.50 

eps. 

Fig. 2. The ellipticity statistics at 2.5 arcmin scale for 
as/0.6 = 2. The statistics for the field points is given in the 
top panels and for the extrema points in the bottom ones. The 
thin lines correspond to the primordial temperature map, the 
thick lines to the lensed map. The error bars correspond to the 
cosmic variance computed for a 9x9 square degree size map. 
The cosmic variance has been estimated with 20 different re- 
alizations of both the primordial temperature maps and the 
displacement fields. The continuous lines give the theoretical 
predictions (34) and (44) with a K given by (48). The dashed 
lines correspond to the best fit of the ellipticity distribution 
using formulae (34) or (44) where a K is assumed to be a free 
parameter. 
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where 5r(i, j) is the local temperature fluctuations at the 
grid point (i,j). The Figs 2-3 present the results of the 
distribution function for various filtering scales, and var- 
ious normalization of the lens effects (expressed in terms 
of as). We can see that for a very small filtering scale, 2.5 
arcmin, the numerical results are in good agreement with 
the theoretical predictions for both the field points and 
the extrema. 

At larger angular, the filtering effect starts to play a 
major role and significantly reduces the lens effects. Thus 
at 5 arcmin angular scale, the amplitude of the effect seems 
indeed to be reduced by a factor of about 2 compared 
to what would be expected when the smoothing effects 
are neglected (see table 1). This discrepancy is shown in 
figure 3 where the dashed lines correspond to the best 
fit of the lens effect when a K is a free parameter. The 
quality of the fit shows however that the filtering effect 
changes the amplitude of the lens effects, but not much 
the e dependence of the corrective terms. 

4-2. The Cosmic Variance 

Because of the use of compensated filters, it is fair to as- 
sume that the a whole sky survey would simply correspond 
to a given number of independent 9x9 square degree maps. 
This should be particularly accurate for the local statisti- 
cal indicators considered here. The cosmic variance then 
scales as the square root of the number of such maps. The 
signal to noise ratio for the lens effect is thus expected to 
obey the scaling property, 



Fig. 3. Same as figure 2 
crs/0.6 = 4 



with a 5 arcmin filtering scale, Signal 



oc 



Noise \/Sky coverage 



(52) 



The maps are filtered with a compensated filter built 
as the difference of two Gaussian window functions with 
one filtering radius being twice the other, 

W 8o (l) =cx P H 2 o 2 /2) -cxp(-2Z 2 2 ). (49) 

This make the maps rather weakly sensitive to the large 
angular scale correlation. The correlation length of the 
local quantities is thus expected to be rather small, of the 
order of 2 8q. 

The local curvature is calculated from finite differ- 
ences. To be more precise the scheme which is used is 
the following, 

ST,u(i,j)= ^ S T {i-l,j + k)- 

fe=-i,i 

2S T (i,j + k) + 8 T {i + l,j + k); (50) 
*r,i2(M') = \ [*r(* 6 T {i - l,j + 1)- 

6 T (i + l,j-l) + 6 T (i + l,j + l)]. (51) 



If one is able to get a the sky coverage of 50% as an- 
ticipated for the Planck surveyor then it means that the 
cosmic variance should be, for the whole survey, about 16 
times less important than for a 9 x 9 square degree map. 
The signal to noise ratio is expected to be the same for a 
"whole" sky survey with as — 0.6 or for a 9 x 9 square de- 
gree map with as = 4 x 0.6. The latter case then provides 
us with a realistic estimation of the Cosmic Variance for 
the future Planck mission. 

In table 1 , I summarize the results obtained for the de- 
termination of a K with a best fit and the error associated 
with this determination. This error has been determined 
as the r.m.s. of the determined cr 2 among the 20 realiza- 
tions at my disposal. One can see that the signal ratio 
grows indeed with as and scales like a\ for a given filter- 
ing scale. The filtering effects however reduce the signal 
by a factor 3 from what would be expected with a direct 
determination of a K . The noise is independent on as but 
depends slightly on the filtering scale. The signal to noise 
ratio expected for the Planck survey at 5' resolution scale 
is thus about 3 if as = 0.6. 
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Table 1. Values of a K from direct determination with the displacement maps (eq. 48) and from best fit parameterizations. 
Errors give the cosmic variance for maps of 9x9 square degree size. 



a 2 K X 10 3 

as ang. scale direct det. 



0.6 
1.2 
1.2 
2.4 
2.4 



o'i x 10 3 , best fit, 
field points, 
no lens effects 



al x 10 3 , best fit, 
field points, 
with lens effects 



al x 10 3 , best fit, 
extrema, 
no lens effects 



2.5' 
2.5' 

5' 
5' 

+ noise 



2.3 ±0.1 
9.2 ±0.3 
5.1 ±0.3 
20. ± 1. 
20. ± 1. 



0.2 ± 1.2 
0.2 ± 1.2 
0.4 ± 1.6 
0.4 ± 1.6 
0.4 ± 1.5 



2.5 ± 1.4 
9.4 ± 1.5 
1.8 ± 1.6 

5.6 ± 1.6 
5.6 ± 1.6 



0.2 ± 1.8 
1.2 ± 1.8 
1.0 ± 2.1 
1.0 ± 2.1 
1.2 ± 2.0 



al x 10 3 , best fit, 

extrema, 
with lens effects 



4.0 ±1.7 
10.3 ±1.8 
2.5 ± 2.1 
6.5 ± 2.2 
6.5 ±2.1 



4.3. The effect of Noise 



0.0 0.2 0.4 0.6 0.8 1.0 

eps. 



Fig. 4. Same as figure 2 with a 5 arcmin filtering scale, 
as/0.6 — 4 and including noise 

Fig. 4 and the last line of the table show the effects of 
noise of the measured map. Its effect is found to be very 
weak. However the noise here has been assumed to be a 
pure white noise, and thus, for the ellipticity statistics, 
it reproduces the Gaussian case. This is certainly not a 
realistic assumption. The striping introduces indeed non- 
local correlations that may affect significantly statistics 
that are related to local topological quantities. The study 
of the striping effect is very dependent on the experiments 
and is beyond the scope of this paper. 

5. Discussion 

5.1. The dependence on the cosmological parameters 

The dependence on the cosmological parameters is ex- 
plored here in a quite rough way. The temperature power 
spectrum, and its relation to the density power spectrum, 
has indeed a rather complicated dependence on all the cos- 
mological parameters. So in this simple study I will assume 
that the shapes of C; and P{k) remain the same and sim- 
ply discuss the dependence of the amplitude of o K on the 




0.0 0.2 0.4 0.6 0.8 1.0 

Omega 



Fig. 5. The ratio a 2 K (Q, A)^ 1 /a 2 K (n = 1, A = 0) as a function 
of and A for a power law spectrum with n = — 1.5. 

cosmological parameters due to the density-convergence 
relationship. The ratio o K jo% is obviously a quantity that 
one might want to consider but it should be compared to 
other constraints coming from large-scale structure forma- 
tion. The number density of clusters can provide us with a 
constraint which is in principle free of bias contamination. 
It constrains the amplitude of the density fluctuations at 
roughly the 8 h^ 1 Mpc scale for an Einstein-de Sitter Uni- 
verse (Oukbir et al 1997), which corresponds roughly to 
the scale of interest for the lens CMB effects. For open 
universes the dependence cannot be given simply in terms 
of as because linear mass scale of the galaxy clusters is 
shifted to a larger scale. A dependence arises then with 
the slope of the power spectrum (Oukbir & Blanchard 
1997). Ignoring these subtleties, and following Eke, Cole 
& Frenk (1996) I will simply assume that the actual con- 
straint coming from the observed number density of rich 
clusters can be written, 

a$ Q - 5 w 0.6 ±0.1. (53) 

It means that the larger fl the smaller cr 8 . To estimate 
the variation of a K with roughly a fixed number density of 
clusters, one should then compute er;?(f2,A) ^l~ 1 /a1{Vl = 
1,A = 0). The result is given in Fig. 5. One can see, as 
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expected, that the magnitude of the signal is very weakly 
ft dependent if A = but significantly grows with A. It 
implies that the signal to noise ratio would be a factor 
about 2 larger for a A-CDM. 

5. 2. On the interest of lens effect detection. 

The observational quantities that have been investigated 
here have been designed to be sensitive to the non-linear 
couplings induced by lens effects. It is of course possible 
to consider other quantities, such as the Minkowski func- 
tional, whose general properties for a Gaussian field have 
been recently investigated in detail (Winitzki & Kosowski 
1997 and Schmalzing & Gorski 1997). The lens effects 
might indeed significantly affect the Minkowski functional 
for 2D maps. For instance the shapes of high threshold 
peaks are shown to be more elongated, which implies that 
the averaged circumference near the top of the peak should 
increase. However, a complete theoretical investigation of 
this effect is quite difficult because this is a nonlocal in- 
dicator. And in general, the most efficient non-Gaussian 
indicators probably depend on the processes one wants to 
detect. The number of hot and cold spots seems to be a 
good way to detect topological defects. Is is obviously not 
the case for the lenses. 

In a previous paper the four-point function was con- 
sidered. Actually the local curvature can be viewed as an 
other way to have access to the same information. Here the 
four point function is simply averaged in a way that avoid 
too much cancellation. And other possible quantity is the 
collapsed four-point function, that identifies with the four 
order cumulant of the local temperature. This quantity 
however tends to cancel the contribution of the various 
terms. This seems not to be the case for the local curva- 
ture. It does not mean however that we cannot do better 
for detecting the lens effects. For instance the two-point 
correlation function of the local ellipticity (with orienta- 
tion) might be a good indicator. However one should have 
in mind that the number of CMB structures per lens is 
rather small. It is thus not obvious that it can improve 
the situation. 

Note that the results obtained for the lens effects, tem- 
perature four-point function or ellipticity statistics, are 
all sensitive to the lens two-point correlation function. It 
comes from the fact that the only quantity related to the 
lenses which is potentially available from CMB maps is its 
power-spectrum. The non-Gaussian properties of the lens 
population are for instance not accessible. 

If it is actually possible to detect a lens effect, the 
smallness of the signal to noise ratio with which it can 
be determined indicates that it will be pointless to use 
this information as a mean to constrain the cosmological 
parameters. But, it would be extremely interesting to be 
able to check that the amount of lens effects computed 
from models favored by the CMB power spectrum arc in 
agreement with its detection, or its non-detection. 



Note however that other secondary effects might as 
well induce non-Gaussian properties. The Rees-Sciama ef- 
fect in particular is due to the evolving non-linear po- 
tentials and has a source term which is intrinsically non- 
Gaussian. A peculiar case of this effect is the moving lens 
effect (Birkinshaw & Gull 1983): potential wells that move 
perpendicularly to the line of sight induce temperature 
fluctuations. This effect is quadratic with the cosmic fields 
(proportional to the local velocity times the gradient of the 
potential) and is thus intrinsically non-Gaussian. These 
cases are probably worth investigating. I expect however 
that they give a smaller non-Gaussian signal since the 
static lens effect is the only mechanism that couples the 
primary anisotropies to the line of sight potentials. 
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